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1. Introduction 

Inspiralling and merging neutron star (NS) and/or black hole (BH) binaries are promising sources for 
continuous gravitational waves (GW). Ground-based laser interferometers as e.g. LIGO, VIRGO, and GEO 
are already searching for those astrophysical sources [Tj. For a successful search with the help of matched 
filtering of the emitted GW signals, one needs a detailed knowledge of the orbital dynamics. Spin effects of 
higher order were discussed in [21 El HI El IS] for the inspiral of compact binaries were the orbits were assumed 
to be quasi-circular. A recent publication [7] gave a numerical insight into the evolution of binary systems 
having spins that are parallel to the orbital angular momentum and evolving in quasi-circular orbits. 

Because there are many physical degrees of freedom involved, it is computationally desirable to have 
an analytical description, especially for interferometers working in the early inspiral phase, where numerical 
relativity currently fails to produce hundreds of orbital cycles. For non-spinning compact binaries, the 
post-Newtonian (PN) expansion in the near-zone has been carried out through 3.5PN order [5] and 3.5PN 
accurate inspiral templates have been established for circular orbits UOj . For numerical performances of 
these templates see [TTJ [12]. Observations lead to the assumption that many astrophysical objects carry 
a non-negligible spin, such that the effect of spin angular momentum cannot be ignored for detailed data 
analysis. The problem of spins in General Relativity (GR) was first discussed in |13lfHl[T5] and considerable 
further developments were made in the 1970s [TH1 El Q~5J Ql|] , and in recent years as well. Apostolatos [20] 
showed in his analysis of simple precession for "circular" orbits and spinning self-gravitating sources that the 
form of the GW signal is affected. The amount of the energy radiated by the binary system with spin has 
been determined by [21| . 

Therefore, we want to implement recent breakthroughs in dynamical relativity of spinning compact 
binaries into a useful prescription for data analysis applications. Our aim is to connect the following items. 

(i) The "standard" procedure to compute the evolution of eccentric orbits from the Hamilton equations of 
motion (EOM). For eccentric orbits, but neglecting spin effects, Damour and Deruelle [22] presented a 
phasing at 1PN employing conchoidal transformations to reduce the structure for the radial motion. 
Later publications [231 HH [25] used Hamilton EOM instead of Lagrange ones and employed a more 
general scheme for a solution to conservative 3PN dynamics without spin. 

(ii) The 2PN point particle (PP), next-to-leading order spin-orbit (NLO-SO) and leading-order spin(l)- 
spin(2) (LO-S1S2) and spin(l)-spin(l) (LO-S 2 ) interaction contributions. 

As a starting point, we assume (anti-) aligned spin and orbital angular momentum vectors for an estimation 
of the effects. It is interesting to analyse this matter system configuration because numerical results of a 
recent publication indicate that maximum equal-spins aligned with the orbital angular momentum lead to 
observable volume of up to ~ 30 times larger than the corresponding binaries with the spins anti-aligned 
to the orbital angular momentum |26| . From Figure 10 in |26) . one can also find an observable volume of 
those binaries up to ~ 8 times larger compared to non-spinning binaries. These authors conclude that those 
systems are among the most efficient GW sources in the universe. In another recent publication |27) it can 
be found that in gas-rich environments the spins of two black holes can align with the larger scale accretion 
disc on a timescale that is short as 1% of the accretion time. Due to the model of those authors, having two 
black holes interacting independently with an accretion disc, their spins tend to be aligned with each other 
and with the orbital angular momentum more or less depending on the model parameters. 

We work only with the conservative Hamiltonian for the time being, and restrict our attention to terms 
up to 2.5PN order overall, assuming maximally spinning holes. This means neglecting both the well-known 
3PN PP contributions, and the NLO-SiS^ [25]) as well as the NLO-S* 2 contributions, which have recently been 
derived for general compact binaries [29]. This latter publication came out at a late stage in our calculations, 
but it should be a straightforward task to include these terms in a future publication. 

If the objects are slowly rotating, the considered leading-order spin-squared contributions are shifted to 
3PN order and, for consistency, the 3PN PP Hamiltonian has to be included. The 3PN PP contributions to 
the orbital elements are available in the literature |25| and simply have to be added to what we are going to 
present in this paper. Anyway, this work is consistently worked out to all terms up to 2.5PN, having maximal 
rotation or not and will list all results in the spins which are counted of OPN order. 

The paper is organised as follows: Section [2] summarises and discusses the Hamiltonian terms we want 
to include in our prescription. Section [3] investigates the conservation of initial spin and orbital angular 
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momentum alignment conditions. In Section 01 we briefly summarise the Keplerian parameterisation for 
Newtonian orbital dynamics and outline the generalisation to higher PN accurate dynamics. The solution of 
the Hamilton EOM is the subject of Section [5] Section [6] summarises all important results of our procedure. 
In Section [7J we give some formulae for the polarisations of the gravitational waves which are emitted by the 
system. Calculations were mostly done with Mathematica and xTensor (see [301 131) and references therein). 



2. Spin and orbital dynamics 

In the following sections, the dynamics of spinning compact binaries is investigated, where the SO 
contributions are restricted to NLO and the S1S2 and S 2 to LO. The PP contributions are cut-off after 
the 2PN terms. The Hamiltonian associated therewith reads 

H(xi,X 2 ,pl,p2,Sl,S2) = Hpp + Hpp N -(- Hpp N + HgQ + HgQ + HgP + Hg lS2 . (1) 

These are sufficient for maximally rotating black holes up to and including 2.5PN. The variables p a and x a 
are the linear canonical momentum and position vectors, respectively. They commute with the spin vectors 
S a , where "a" denotes the particle label, a = 1,2. H pp is the conservative point-particle ADM Hamiltonian 
known up to 3PN, see, e.g., [32] and [33] . The LO spin dependent contributions are well-known, see, e.g., 
[I3CEH1EI]. H SO° was recently found in [3J3 [35] and H^ L S ° in [3J[2S] (the latter was confirmed in [37]). The 
leading-order S 2 and S 2 . Hamiltonians were derived in [18] and [38]. Measuring the GW signal, determinantion 
of constraints of the equation of state of both extended bodies is possible in principle. Hamiltonians of cubic 
and higher order in spin are given in (39] [40], and higher PN orders linear in spin are tackled in |41[ [42] . 

The four-dimensional model behind the Hamiltonians linear in the single spin variables is given by the 
Mathisson-Papapetrou equations |14[ [T3] 

2PM , (2) 



dr 

Drf 1 



dr 2 pPa 

and the Tulczyjew stress-energy tensor density 



^p0a U a^a a ! (3) 



(4) 

which can be used as the source of the gravitational field in the Einstein equations (see [45] and references 
therein for spin-squared corrections in the stress-energy tensor). Here, the four-dimensional coordinate of the 
a-th object is denoted by x% and p api is the linear momentum, the 4- velocity, normalised as u^u^ = — 1, 
t the proper time parameter, S£ v the spin tensor, ";" denotes the four-dimensional covariant derivative, and 
<5(4)a — S(x — x a ) with normalisation J d^x 5^ a = 1. 'Bf i p p a is the four-dimensional Riemann tensor and D 
the four-dimensional covariant derivative. In order to close the system of equations, one has to impose a spin 
supplementary condition (SSC), which is most conveniently taken to be 

S% V p av = 0. (5) 

To linear order in spin, p afl = m a u a ^, where m a is the mass parameter of the a-th object. Notice that the 
matter variables appearing in the Mathisson-Papapetrou equations and the stress-energy tensor are related 
to the canonical variables appearing in the Hamiltonians by rather complicated redefinitions. 

We are going to the centre of mass (COM) frame, where the total linear momentum vector is zero, i.e. 
p>2 = — pi = —p. The Hamiltonians taken into account depend on ici and X2 only in the combinations 
xi — x 2 , so they can be re-expressed in terms of 1112 = — 1121 = x/f = ~x.fr, the normalised direction from 
particle 1 to 2, and f = |xi — X2I with x = xi — X2. 

We will make use of the following scalings to convert quantities with hat to dimensionless ones, 

H = , (6) 

IJ,C Z 

Gto n 



p = p {jicy 1 , (8) 
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S a = S a ^ — ^(m a c)| , (9) 

h = rni2 x p. (10) 

Here, m = mi + ma denotes the total mass and /i = mi 7712/771 is the reduced mass. The speed of light is 
denoted by c and G is Newton's gravitational constant. 

Explicitly, the contributions to the rescaled version of Equation ([1]) read 

~2 i 



HIM f 
PP - — 



ttIPN 
Hpp 



r(3r?-l)( P 2 ) 



2\ z 



l r 



(3 + r/)(p 2 )+r/(ni 2 ■ p) 



1 

2^ 



H 



2PN 
PP 



— (1-577 + 577 2 ) (p 2 )' +- 



(5 - 20n - 3r? 2 ) (p 2 )' 



2/7 2 (nia ■ pY (p 2 ) - 3r/ 2 (n 12 • p) 4 



1 

2 L 



(5 + 8??) (p 2 ) +3 ?? (ni 2 -p)' 



H 



LO 
SO 



+ v / T^)(h-Si)+ (i-|-yr3^) (h . S2 )} 



H 



NLO 
SO 



e 4 <5 



16r 4 



(h-Si 



12?yr 1 - 77 + Vl - 4?7 (ni2 ■ p) 



J] r (9 - 677 + 19-y/l -477) (p 2 ) 
16 ((ry + 3)Vl - 4jj + 3] 



(h-S 2 ) 



12?/r ^-1 + 77 + v/1 - 4 ?? J (ma ■ p) 
+ J] r (-9 + 677 + 19 v/1 - 477) (p 2 ) 
-16 (fa + 3)^1^- 3) |, 



LO 

S1S2 

LO 



2 r2 

e a 



e 2 <5 2 



ft -j{3(ni 2 -Si)(nia-Sa)-(Si-S 2 )} , 
2r 3 



Ai 



A- 



1 + 2n - ^T^n) (3 (ma ■ Si) 2 - (Si ■ Si 
1 + 27/ + v/I - 4^) (3 (ma • S 2 ) 2 - (S 2 • S 2 )) | . 



(11) 
(12) 



(13) 
(14) 



(15) 
(16) 

(17) 



We introduced dimensionless "book-keeping" parameters e to count the formal 1 / c order and 8 to count the 
spin order (linear or quadratic). Evaluating all given quantities, those have to be given the numerical value 1. 
The parameters a so , ct SlS2 , a s i distinguish the spin-orbit, spin(l)— spin(2) and the spin-squared contributions 
and can have values 1 or 0, depending on whether the reader likes to incorporate the associated interactions. 

The spins are denoted by Si for object 1 and S2 for object 2. Notice that the S 2 and S 2 Hamiltonians 
depend on constants Ai and A 2 , respectively, which parametrise the quadrupole deformation of the objects 
1 and 2 due to the spin and take different values for, e.g., black holes and neutron stars. For black holes, 
X a = — i and for neutron stars, A a can take continuous values from the interval [—2, —4] [3H 06J Q. 

The parallelism condition tells us to set the spins to S a = XcM/h, where —1 < Xa < 1- During our 
calculations, we insert the condition of maximal rotation (S a ~ e) to cut off every quantity after 2.5PN, but 
list our results in formal orders S a ~ e° (for the formal counting, see, e.g., [3D] and also Appendix A of |42j). 



I Note that the definition of the A a depends on the definition of the spin Hamiltonian and, thus, can be arbitrarily normalised. 
We consistently use the notation mentioned above. 
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However, for S a ~ e 2 , many spin contributions are of the order 0(e 6 ), i.e. 3PN which is beyond our present 
2PN PP dynamics. The reader may insert either S a ~ e (maximal rotation) or S a ~ e 2 (slow rotation). The 
next step is to evaluate the EOM due to these Hamiltonians and to find a parametric solution. As stated, 
we will restrict ourselves to parallel or anti-parallel angular momenta and will, finally, only have to take care 
of the motion in the orbital plane. 



3. Conservation of parallelism of h and Si and S2 

The motion of binaries with arbitrarily-oriented spins is, in general, chaotic as soon as the spin-spin interaction 
is included [37125]. For special configurations, despite this, it is possible to integrate the EOM analytically, 
which particularly is the case for aligned spins and orbital angular momentum. 

The time derivatives of the spins S a and the total angular momentum J are governed by the Poisson 
brackets with the total Hamiltonian, given by 

[Si,H] = <5 2 e 2 ja SlS2 -J (3 (S 2 ■ n 12 ) (n 12 x Si) + (Si x S 2 )) 
(n i2 x Si) 



(Si • nia) 3Ai (27? - 1 - y/l - 4t?) J 

+ *{^<!^>(-| + v^+i) 

(h * Sl> ^(p-n 12 ) 2 ,,(l->, + % /T^) 



+ i(P> (9 - 677 + 19^1-477 



16 

(hxSi) 



(3 + (77 + 3)yi _ 4^) 



[S a ,H] = [Si,H](lo2), (19) 
[J,H] = [h,H] + [Si,H] + [S 2 ,H]=0. (20) 

Furthermore, the magnitudes of the spins are conserved, because the spins commute with the linear 
momentum and the position vector and fulfill the canonical angular momentum algebra. Note that the 
operation (lf>2) switches the label indices of the individual particles and goes along with ni 2 o n 2 i = — ni 2 . 
Equation f|20[) is not displayed completely here like Equation (fT8j) . If we assume parallel spins and orbital 
angular momentum at t = 0, all the above Poisson brackets vanish exactly. Anyway, this is insufficient to 
conclude the conservation of parallelism of h and the spins for all times t > since 

Si(t) = Si|t= + [Si,H] \ t=0 t+ \ [[Si,H] ,H] \ t=0 t 2 + ... 



OO 1 



n 

n=0 



H]Jt=ot n , (21) 



where 



[Si,H]„ = [[Si.H^^H] , 

[Si,H] =Si. (22) 

Because the system of variables Si and S 2 has to be completed with r and p to characterise the matter 
system, one has to give clear information about the full system of EOM. It is important that, even with 
vanishing Poisson brackets of H with Si and S 2 , r and p do change due to the orbital revolution. Thus, 
one has to clarify if this non-stationary subsystem of the EOM is able to cause violation of the parallelism 
condition during time evolution. From the stability theory of autonomous ordinary differential equations it 
is well known that there is a fixed point if all time derivatives of the system vanish. In the case of a system 
starting at such a fixed point at t — it will not be able to evolve away from this point. The discussion of 
these issues is the main point in the following two subsections. 
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3.1. Discussion via conservation of constraints 

One way to show the non- violation of the initial constraint of h || Si, S2 due to the motion of the binary is 
to argue via the time derivatives of the constraints. These should be written as a linear combination of the 
constraints themselves. Let 

C a (x,p,S) = 0, (23) 

be the initial constraints of the system. Dirac [491 p. 36] argued: If one can write 

C a =J2 D ab( x ,P,S)C b , (24) 

b 

for the time derivatives of the constraints, the constraints are conserved. That is due to the fact that every 
time derivative of Equation (|24[) generates only new time derivatives of the constraints on the one hand, 
which can be expressed as a linear combination of constraints, or time derivatives of the quantities appearing 
in D a b times the constraints on the other. 
In our case the constraints read 

Sx-^±ih=Si-xih = 0, (25) 

S 2 -!^h=S 2 -x 2 h = 0, (26) 
h 

with xi aid Xi denoting the ratios of the spin lengths and the orbital angular momentum. In general, the 
quantities Xa have non-vanishing time derivatives, 

hh 



i -t-- , ¥M-^' ™ 

Due to the conservation of the total angular momentum J, the derivatives of ([23)) and (f2"6")l can be expressed 
via 

d _ . , dS a _ ( h ' h ),_ . dh 

-(S Q - Xa h) =_ +Xa ^^h- Xa - 

dS a / h h\ dh 
dS Q - / h h\ ^ dS b 

b 

where the tensor product l/h 2 h<£> h is the projector onto the h direction. Note that the constraint equations 
only depend on the spin derivatives in linear manner. Hence, it is sufficient to analyse the structure of 
Equation ([T5|. 

^ =D 1 (S r n 12 )(n 12 xS 1 ) + D 2 (S 1 xS 2 ) 

at 

+ D 3 (Si • ma) (n 12 x Si) + D 4 (hx Si) . (29) 

The coefficients Dk (k = 1, ... ,4) are all scalar functions of the linear momentum p, the separation r and 
other intrinsic quantities. We are allowed to add vanishing terms to Equation (|29p . namely 

^ =D 1 [(S a • ni 2 ) -£2 (h- n 12 )] (n 12 x Si) + D 2 (S x x S 2 ) 

+ D 3 [(Si • ma) - xi (h • ma)] (ma x Si) + D 4 (h x S x ) . (30) 
As well, we can add a term to the D2 coefficient and subtract it at the end, getting 
dS 

— - = Di [(S 2 • nia) - Xi (h • ni 2 )] (ni 2 x Si) 

dr 

+ D 2 [(Si x S 2 )-xa(Si x h)] 

+ D 3 [(Si • ma) - Xi (h • ma)] (n i2 x Si) 

+ (D 4 -D 2 x 2 )(hx Si) . (31) 
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Finally, we can insert a vanishing term into the modified last one: 
dS 

— - = L>i [(S 2 • ma) -X2 (h • n i2 )] (ni 2 x Si) 

dr 

+ D 2 [(Si x S 2 )-x 2 (Si x h)] 

+ D 3 [(Si ■ m 2 ) - xi (h ■ m 2 )] (ma x Si) 

+ (D 4 -D 2X2 )[(hx Si)-xi(hxh)] . (32) 

We still need to compute the time derivative of S 2 to obtain the full derivative of the constraints. Therefore, 
let Ek be the scalar coefficients in dS 2 /di (equivalent to the Dj~ in ([521 ). Using this, we can rewrite it as 

^ = Ei [(Si • n i2 ) - xx (h ■ ma)] (ni2 x S 2 ) 

dr 

+ E 2 [(S 1 x S 2 )-xi (hx S 2 )] 

+ E 3 [(S 2 • m 2 ) - X2 (h • ma)] (n i2 x s 2 ) 

+ {E A + E 2X i)[(hxS 2 )-X2(hxh)} . (33) 
Thus, the complete time derivative of e.g. the Si constraint (|28|) is given by 

S(Si-Xih) =((l + fe )l-Xi^^ 

£>i [(S 2 • ma) - x 2 (h • n 12 )] (n 12 x Si) 

f D 2 [(Si x S 2 ) - X 2 (Si x h)] 
f D 3 [(Si • ma) - Xx (h • ma)] (m 2 x Si) 



+ (fl 4 -fl 2 X 2 )[(hx Si) -xi (hxh)] 
/ h h 

£?i [(Si • ma) - Xx (h • ma)] (ni 2 x S 2 ) 

+ S 2 [(Si x S 2 ) - xi (h x S 2 )] 
+ £ 3 [(S 2 • ma) - X2 (h ■ ni 2 )] (m 2 x S 2 ) 

+ (E 4 + E 2 xi) [(h x S 2 ) - X2 (h x h)]l . (34) 



fn each of the summands of Ek and 1?^ in the above equation, one can factor out the constraints linearly. 
Thus, they do vanish if the constraints are inserted. 

3.2. Discussion via symmetry arguments 

To underline the results of subsection 13 . 1 1 we want to show that the multi Poisson brackets ([22[) all vanish if 
we demand the parallelism of h and the spins, as a complement to the constraint evolution analysis. During 
the calculation of the expressions, we truncated the terms to quadratic order in spin and to 2.5PN order, 
counting the spin maximally rotating. There is a finite set of terms which are axial vectors and linear or 
quadratic in spin. Additionally, the spin has to appear in a vector product. The reason is that the spins 
commute with the PP Hamiltonians, and Poisson brackets of spins with spin Hamiltonians will give cross 
products of spins with angular momentum h or Si, S 2 , respectively. In the Hamiltonians, there are only 
scalar products of spins with other vectors or with the spins themselves, so that the e.yfc are still remaining 
after evaluation of the Poisson brackets. The products with the correct symmetry and linear and quadratic 
in spin are of the form 

A • ^ and^-^\S 

Ax A A- A axp a-p 

P X P P P 
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where A stands for axial vectors, P for polar vectors, S for scalars and PS for pseudo scalars. The vector 
A can be h, Si, S2 which are axial vectors and P can be ni2, p which are polar vectors. Now we enumerate 
all spin products that emerge when we compute the first multi Poisson brackets, omitting scalar factors like 
e.g. functions of r and 77: 

[Si, H] : (Si x n 12 ) (Si ■ n 12 ) ; (Si X n 12 ) (S 2 • n 12 ) ; 

(SixS 2 );(hxSi) , (35) 
[[Si,H],H] : (Si x m 2 )(p-Si);(Si x n 12 ) (p • S 2 ) ; 

(Si x p) (ma • Si) ; (Si x p) (n 12 • S 2 ) , (36) 
[Si,H] 3 : (SiXp)(p-Si);(SiX P )(p-S 2 ); 

(Si x p) (ma • Si) ; (Si x p) (n 12 • S 2 ) . (37) 

Additionally, terms coming from Poisson brackets of S 2 with H appear and can be computed from above 
with the operation (1 <-> 2). Evaluation of higher multi Poisson brackets will not create new terms, but only 
increase the number of already known factors in the products. These terms vanish identically if we consider 
the parallelism between h and the spins because of (h • m 2 ) = and (h • p) = 0, by construction. Due to the 
vanishing Poisson bracket of J = h + Si + S 2 with the Hamiltonian H (this is true even without demanding 
the parallelism), the disappearance of all multi Poisson brackets of Si and S 2 with H turns out to be sufficient 
to conclude the constancy of h if the motion starts with h j| Si, S 2 . 



4. Kepler Parameterisation 

In Newtonian dynamics the Keplerian parameterisation of a compact binary is a well-known tool for celestial 
mechanics, see e.g. [50]. After going to spherical coordinates in the COM, (r,8,(f>) with the associated 
orthonormal vectors (e r ,ee,e<^) and restricting to the 9 = 7r/2 plane, the Keplerian parameterisation has 
the following form: 

r = a (1 — e cosu) , (38) 
<j>-fo=v, (39) 



2 arctan 



1 + e u 
tan — 

1 - e 2 



(40) 



Here, a is the semimajor axis, e is the numerical eccentricity, u and v are eccentric and true anomaly, 
respectively. The time dependency of r and <p is given by the Kepler equation, 

I = n (t — to) = u — e sin u , (41) 

where I is the mean anomaly and n the so-called mean motion, defined as n = =p- with P as the orbital 
period [51j . In these formulae to and </>o are some initial instant and the associated initial phase. In terms 
of the conserved quantities E, which is the scaled energy (see Equation ©) and numerically identical to H, 
and the orbital angular momentum h, the orbital elements e, a and n satisfy 

a =wv (42) 

e 2 = 1 - 2h 2 \E\ , (43) 
n = (2|£|) 3 / 2 . (44) 

For higher PN accurate EOM it is possible to get a solution in a perturbative way, having the inverse speed 
of light as the perturbation parameter. The 1PN accurate Keplerian like (from now on we refer to quasi- 
Keplerian) parameterisation was first found in |22] and extended for non-spinning compact binaries in |23| 125) 
to 2PN and finally 3PN accuracy. 

In the recent past a number of efforts has been undertaken to obtain a solution to the problem of spinning 
compact binaries via calculating the EOM for spin-related angular variables in harmonic gauge. For circular 
orbits, including radiation reaction (RR), the authors of |52| evaluated several contributions to the frequency 
evolution and the number of accumulated GW cycles up to 2PN, such as from the spin, mass quadrupole 
and the magnetic dipole moment parts. The gravitational wave form amplitudes as functions of separations 
and velocities up to and including 1.5PN PP and 1.5PN SO corrections are given in [S3], discussed for the 
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extreme mass ratio limit in the Lense-Thirring approximation and later in |54| and |55| for comparable mass 
binaries. Recently, in |56| a set of independent variables and their EOM, characterising the angular momenta, 
has been provided. 

For circular orbits with arbitrary spin orientations and leading-order spin-orbit interactions, the spin 
and orbital solutions for slightly differing masses were given in |57] . Including LO contributions of S 2 , 
S1S2 and SO as well as the Newtonian and 1PN contributions to the EOM, a certain time-averaged orbital 
parameterisation was found in |58| , for a time scale where the spin orientations are almost constant, but 
arbitrary and the radial motion has been determined. Symbolically, those solutions suggest the following 
form for the quasi-Keplerian parametrisation including spin interactions: 

r = a r (1 — e r cos u) , (45) 

n(t — to) = u — e t sinii + F v - U (v — u) + T v sinw + Ti v sin2v + T§ v sin3w + . . . , (46) 



2tt 



- (j) ) = v + Qiv sin2ii + G 3v sin3v + Gav sin4u + Q 5v sin5w + . . . , (47) 



2 arctan 




(48) 



The coefficients J 7 ..., C*... are PN functions of E, h and r\. At the end of the calculation for binary dynamics 
with spin, they will obviously include spin dependencies as well. 

In case RR is included, the orbital elements are not longer to be regarded as constants. Damour, 
Gopakumar and Iyer published equations of motion for these elements for the case that the RR is a small 
effect and the time derivatives due to RR will not contribute to the GW expressions explicitly [59l [60] . Spin 
effects of RR in eccentric orbits were discussed in (6T1 and references therein. 



5. The quasi-Keplerian parameterisation for aligned spinning compact binaries 

Having proven constancy in time of the directions of angular momenta, we can adopt the choice of spherical 
coordinates with h || eg (in the 9 = it/2 plane) and the basis (ni2 = e r , e^). Hamilton's equations of motion 
dictate 

H 

r = n i2 • r = n i2 ■ , (49) 
op 

<9H . , 

r<p = es ■ r = es ■ —— , (50) 
op 

with r = dr/dt and <p = d(p/dt, as usual. The next standard step is to introduce s = 1/r, such that 
r = —s/s 2 . Using Equations (|4"9")l and ([5U|) . we obtain a relation for r 2 and thus s 2 and another one for 
4>/ s = d<p/ds, where the polynomial of s 2 is of third degree in s. To obtain a PN accurate parameterisation, 
we first concentrate on the radial part and search for the two nonzero roots of s 2 = 0, namely s + and s_. 
The results, to Newtonian order, are 



1 _i + vT^F @+0(j2)i 



a r (l — e r ) h 2 



-IZl^S+Oi*). (52, 
a r (l + e r ) nr 

s_ representing perihelion and s + as the apastron. Next, we factorise s 2 with these roots and obtain the 
following two integrals for the elapsed time t and the total radial period P, 

p = 2 r Wdr 

7 S _ tV(t-s-)( S+ -t)' 
which is a linear combination of integrals of the type 
f s+ T n dr 

I' n = 2 - n % =. (54) 



T 



y/(T-S-)(s+ - T) 



The time elapsed from s to s + , 
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is a linear combination of integrals of the type 



h, = 



Mr 



. (56) 

r V( r - s -)( s + - T ) 

Both integrals /„ and I' n are given in |Appendix A| in terms of s + and s_ for /' and in terms of and 
v for /, respectively. The function V^s) is a third order polynomial in s and the factor 2 follows from the 
fact that from s_ to s + it is only a half revolution. With the help of the quasi-Keplerian parameterisation 

r = a r (1 — e r cos it) , (57) 

where a r and e r are some PN accurate semi-major axis and radial eccentricity, respectively, satisfying 

_ 1 s+ + .s_ 

a r ~ j 

2 s_ s + 



1 



(58) 
(59) 

due to (fSTj) and (|52|) . we obtain a PN accurate expression for a r and e r in terms of several intrinsic 
quantities. With Equation (|55p. we get a temporary expression for the Kepler Equation, as we express 



2 s_ 



n (t - 1 ) 



-(t ~ ^o) m terms of u, and as standard, we introduce an auxiliary variable 



2 arctan 



jl + e r 


u 




tan — 


V 1 - e r 


2 



(60) 



At this stage, we temporarily have 
I = n (t - t ) 

= u + T u sin u + Fy- U {v — u) + T% sin{i , (61) 

with J 7 ., as some PN accurate functions of E, h, rj, X a and \a- These functions are only temporarily needed 
and lengthy, so we will not provide them. 

Let us now move on to the angular part. As for the time variable, we factorise the polynomial of dcf>/ds 
with the two roots s_ and s + and obtain the elapsed phase at s and the total phase $ from s_ to s + . 



$ = 2 



V(s--t)(t 
Bs(t) 



-At. 



idr . 



Using Equation (|62f and Equation (|63| . the elapsed phase scaled by the total phase ^ 
v is computed as 

2ir ~ 

— ((/) — <f) ) = v + Gv sin v + Qiv sm 2w + Gzi, sin 2>v . 
For the following, we change from the auxiliary variable v to the true anomaly 



2 arctan 



'1 



1 



— - tan — 

2 



(62) 

(63) 
in terms of 

(64) 
(65) 



due to Equation (|4"5)) with 

= e r (l + e 2 ci + e 4 c 2 ) , (66) 

differing from the radial eccentricity by some PN corrections C\ and c^. These corrections are fixed in such 
a way that the sinw contribution in =£-{4> — 4>o) vanishes at each order and the lowest formal correction to 
the phase is shifted to 2PN. Therefore, we eliminate u in Equation with the help of (pjSj) and insert the 
result into (|64p . Requiring the sin v term to vanish yields 



smi> 




sin v H — c, 

1 4 \ e 2 



■ sin(2w) 



(67) 
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where c\ and c 2 are at most quadratic functions in b and depend on the intrinsic quantities of the system. 
After determining e^, Equation (|62p takes the form 

— {(j) — 4>a) = v + Q2v sin 2v + Qi v sin 3w . (68) 
With the help of v, we can re-express the temporary Kepler equation (|6"Tj) in the form of 

£ = n (t — to) = u — et sinw + J- V — U (v — u) + J- v sinn . (69) 

Here, e t is the time eccentricity and simply represents the sum of all terms with the factor sin it in £. All the 
orbital quantities will be detailed in the next section. 

6. Summarising the results 

We present all the orbital elements a r , e r , et, e^, n and the functions J 7 ., and Q of the quasi-Keplerian 
parameterisation 

(70) 
(71) 

(72) 
(73) 



r = a,.(l — e r cosw) , 

n(t — to) = 11 — et sinu + J- V - U (v — u) + T v sinu . 



2tt 
If 

v 



4>o) =v + Q 2v sin(2u) + Q 3v sin(3i>) . 
= 2 arctan 



'l + e 



4> 



tan — 

1-e,* 2 



in the following list. For 5 — (remember that d counts the spin order) one recovers the results from, e.g. 

EH. 



2\E\ 



■ 77 — 7 <5 



V 1 - 4? 7 (Xi - X2) + (l - |) (xi + 



/ ( 2 



(xi - Xa) 2 + <V Qv^Wi - A 2 ) + i(l - 2tt) (Aj + A 2 )) 

+ (Xi+X 2 ) 2 (a B 2 Q^l - 4^ (Ai-A 2 ) + 1(1-27?) (A x + A 2 )) - 



+ a s 2 (xi + X2) (xi - X2) ( ^(1 - 277) (Ax - A 2 ) + -y/1 -4t7(Ai + A 2 ) 



+ 6 4 || J B|Q(7 ? 2 + 107 / +1) 



+ 4^(11.7-17) 

1 - 2/i 2 |£;| 



(-6r/ 2 + 1977 - 8) (xi + Xi) + gv 71 - 4 ^/ - 8) ( X i - X2) 
(V - ^ + 8J (X! + xa) + J (32 - 97?) y/l - 47? ( X i - X2) 



+ e 2 <j /7 2 |S| 2 ( 5(7? - 3) + -a so [s^I - 4^ ( X i - X a) + (8 - 4t?) ( X i + X2) 



^2 " 

(xi-X2) 2 (2a SlS2 7? + a s 2 (yi - 477 (Ai - A 2 ) - (2r? - 1) (Aj + A 2 ))) 
(Xi + X2) 2 (y (yi - 47; (Ai - A 2 ) - (2r? - 1) (Ai + A 2 )) - 2a SlS2 v) 



h 2 

+ 



(74) 
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+ a s 2 (xi + xa) (Xi " Xa) (Vl-^ (Ai + A 2 ) - 2(2 V - 1) (Ai - A 2 )) ) 
+ \E\ ^-2(r? - 6) + |a so [4(77 - 2) ( X i + X2) - 8^1-47? (xi - X2) 

- p (Xi - X2) 2 (2a SlS2 r; + a s 2 (y/l - 4ry (Ai - A 2 ) - (2 ?/ - 1) (A x + A 2 )j) 
+ (xi + X2) 2 (« S 2 (Vl - 4r; (Ai - A 2 ) - (2r? - 1) (Ai + A 2 )) - 2a SlS2 r?) 
+ a B 2 (xi + X2) (Xi - X2) (2^/1 - 4t? (Ax + A 2 ) - 2(2 ?? - 1) (Ax - A 2 )) ^ | 



e 4 <^ -V + 55t? - 80 



( 6r? 2 _ 49r? + 80 ) (xi + X2) + (g0 _ ig??) 4^ ( Xl _ X2 ) 



+ 



/i 2 



-2277 + 34 



h 

\E\ 2 [ T? 2 +77 + 26 



(-St/ 2 + 7877 - 64) (xi + X2) + 20 - 477(977 - 32) ( X i - X2) 



71 



(IO77 2 - 70 ? ; + 4) (X! + X2) + 4(1 - 4t7) 3 / 2 (xi - X2) 



e 4 4 



677 - 15 



CH-so , 

h 



2 (t? 2 -877 + 6) (xi + X2) - Ay/1 - 4/7(77 - 3) (Xi - Xi) 



+ -^l-(ll77 2 + 3077 + 555) 



1 - 2h 2 \E\ + e 2 \ \E\ (h 2 \E\(17 - 777) + 4(77 - 1)) 



+ -a so |S| [2(77 - 2) (xi + X2) - 40 - 477 (xi - X2) 



S 2 \E\ 
h 2 



(Xi -X2) 2 (a s 2 (Cn-^j (A1 + A2) - ^0 -4t?(Ai -A 2 )J -a SlS2 ??J 
+ (xi + X2) 2 (^a Sl s 2 r/ + a s 2 ( U - (Ai + A 2 ) - i0-4? ? (Ai - A 2 )^ J 
+ <V (Xi + X2) (Xi - X2) ((277 - 1) (Ai - A 2 ) - 0- 477 (Ax + A 2 )) J 
1 j JjpY-llji + 17 + /j 4 |£f (-16?7 2 + 47t? - 112) + 12\/2h 3 \E\ 3 / 2 (5 - 2n) 



+ 2h 2 \E\ (5?7 2 + 77 + 2) + 6\/2hy/\E\(2r) - 5) 



6 \E[ 
+ h so 2/7 2 



(xi + X2) (-16V2h 3 \E\ 3/2 (t? 2 - 877 + 6) + h 2 \E\ (32t7 2 - 159t7 + 124) 
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+ 8V2 s/\E\h (rf - 877 + 6) - 8?f + 78r/ - 64^ 
+ y^T^ (xi - Xi) (z2^/2h z \E\V 2 (n - 3) + ft 2 |£|(124 - 5 
- lQV2h^J\E\{ri - 3) + 18?7 - 64^ 



•J V — U 



T = 



<i> 

2^ 



Q-iv 



3 [V 



h 

,|£| 3 / 2 
2V2h 



{if - 87? + 6) (xi + X2) - 2^1-4^(7? - 3) (xi - xa) 



VI -2/i 2 |S| -77(77 + 4) 



l + e 



/i 2 



V 1 - 477(77 + 8) (xi - X2) - (1377 - 8) (xi + X2) 

, 5 
3 + -a s 
h 



h 2 



(7/ - 2) (xi + X2) - 2y/l -477 (xi - X2) 
I (Xi - X2) 2 (-2a SlS2 7/ - a s2v /l -4r?(Ai - A 2 ) + ^2(277 - 1) (Ai + A 2 ) 
+ I (Xi + X2) 2 (2a SlS2 77 - a s 2 V /l-4?7(Ai - A 2 ) + ^2(277 - 1) (Ai + A 2 )) 
(Xi + X2) (Xi - Xi) ((2r? - 1) (Ai - A 2 ) - v/1-477 (Ax + A 2 )) J 



3 



h 2 
1 fs 



„ 15 , ^ 
377 - T + -a s 

21 



v> 

15 



2 (t7 2 - 87, + 6) (xi + X2) - Vl - 477(77 - 3) (xi 
3 



4 VI - 477(77 - 8) (xi - X2) - ^ (2t7 2 - 4977 + 56) (xi + X2) 



(2r?-7) 



A (2h 2 \E\-l) f 77 ( 37? -1) 



4h 4 

3-a so 
h 



V 1 - 4 W (Xi - X2) - (?? - l)?? (Xi + X2) 
4 (l-2/7 2 |£;|) 3/2 



8h 4 



3rf 
4 



+ tOso (v - 1)t? (Xi + X2) - V 1 — 47777 (xi - X2) 
= 1 - 2h 2 \E\ + e 2 (h 2 \E\{ V - 15) + 12) 



+ -a so |£| (h 2 \E\ - 1) [8V1 - 477 (xi - X2) - 4(77 - 2) (xi + Xi) 



+ 



S 2 \E\ 
h 2 



(4h 2 \E\-3) 



(X1-X2) 2 (a SlS2 r) + a s 2 Q V 1 - 4? / ( A i - A 2) - - 1) (Ai + A 2 ) 
+ (Xi + X2) 2 (^-a slS2 T7 + a s 2 Q a /i_4t 7 (Ai - A 2 ) - ^(277 - 1) (A x + A 2 ) 
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+ a s 2 ( X i + xa) (xi - Xa) (V 1 ~*n Oi + A a ) - (2? ? - 1) (Ai - A 2 )) 

. -4h 4 \E\ 2 (3r/ 2 - 30?7 + 160) + 4h 2 \E\ (9r/ 2 + 8877 - 16) - 15?? 2 - 232?? + 408) 

oh 

(Xl + Xa) {2h 4 \E\ 2 {80 - 3I77) + 8h 2 \E\ (rf - 3677 + 17) - 3 (rf - 7l V + 64)) 



14 



S_ \E[ 
h aso 2hi 



+ 0^4^ (xx - X2) (-2/i 4 |£| 2 (r; - 80) + 4/i 2 |£|(34 - 15tj) + 33 )? - 192) 



(83) 



For the case that one chooses et instead of the other eccentricities as the intrinsic parameter to be searched 
for in the data analysis investigations, we give the ratios of the other eccentricities with respect to e t : 



-i- =l + ^(8-3^| 



8\E\ 



{i] - 2) (xi + Xa) - 2^/1 - 477 (xi - Xi) 



+ 



P\E\ 
h 2 



1 



1 



(xi - Xa) 2 (a* ( j(2r? - 1) (Ax + A 2 ) - ^ V ^4^(A 1 - A 2 ) j - ^ 
+ (xi + X2) 2 + « s2 Q(2„ - 1) (Ax + Aa) - \ ^T^l (Ax - A 2 ))) 

+ <V (Xi + Xa) (Xi - X2) ( (v - (Ai - A 2 ) - i v 71 - 4? / ( A i + A 2 ) J | 
/i 2 |£| (6ry 2 - 63r) + 56) - 6V2hy/\F\(27j - 5) - llr] + 17 

Vl -4?7 (xi - X a) f /i 2 |-B|(23t7 - 84) + 16v^fc'\/jl|(»? - 3) + 18?/ - 64 



2/i 2 



5 \E\ 
s °hAh 2 



(Xi + X2) (h 2 \E\ (8n 2 - 5577 + 84) + 8^hy/\E\ (n 2 - 8i] + 6) + 8? ? 2 - 7877 + 64) 



(84) 



1 + e 2 \-2( V -4)\E\ 



S\E\ 



h 2 



(rj - 2) (xi + Xa) - 2^1-4^ (xi - Xa) 
1 kxi Xa) 2 (<v ( (f) (Ai + A 2 ) - \^Tr ] (Ax - A 2 )) - a BlBa ^ 
+ (Xi + Xa) 2 ^ S1S2 ?7 + a s 2 Ur] - M (Ax + A 2 ) - i v /l-4?/(A 1 - A 2 )^ ) 
+ Otf (Xi + Xa) (Xi - Xa) ((2ij - 1) (Ax - A 2 ) - ^/l - 4 V (Ax + A 2 )) J 

2h 2 \E\ (lhf - 168?/ + 224) - 48V2hy/\E\(2r) - 5) - 15?/ 2 - 144?/ + 272 

(Xi + Xa) (-3h 2 \E\ (2n 2 - 15 V + 28) - 8y/2hy/\E\ (if - 8n + 6) + 5T? 2 + 135r? - 128) 



+ a s , 



\E\_ 
16h 2 
S_\E[ 
hAh 2 



+ y/l-4r)(xi - Xa) (/i 2 |£|(13?7 - 84) + 16y/2hy/\E\(ri - 3) + 15r; - 128) 



(85) 



7. Gravitational wave forms 



The final form of the GW model consists of two iixgredieixts. The first one is the expression for the far-zone 
radiation field, which will naturally depend on geixeral kinematic quantities describing the binary system. 
The explicit time evolution of these kinematical quaixtities represents the second one. 
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In the first part of this section, we will compute the GW amplitudes h x and h + . We will require the 
associated corrections to the desired order, which are available in the literature in harmonic coordinates, see 
below for references. As well, we will use coordinate transformations from ADM to harmonic coordinates to 
be able to apply the time evolution of the orbital elements in the previous sections, which we have computed 
in ADM coordinates only. 

In the second part, we will provide the latter quantities as implicit functions of time. 



7.1. PN expansion of the gravitational radiation amplitudes 

The transverse-traceless (TT) projection of the radiation field and thus h x and h + , the two polarisations, 
strongly depend on the observer's position relative to the source. To obtain h x and h + , it is necessary to 
give position relations of the orbital plane to the direction where the detector is situated. An observer- 
dependent coordinate system will be helpful to give the time domain waveform expressions in terms of the 
radial separation, orbital angular velocity and the spins. 

We start the calculation by defining the unit line-of-sight-vector N as pointing from the source to the 
observer. Now, let the unit vectors p and q span the plane of the sky for the observer and complete the 
orthonormal basis (p, q, N), 

pxq = N and cyclic. (86) 

Additionally, let us define an invariant reference coordinate system (e x , e y , e z ). Both coordinate systems can 
be coupled by a special orthogonal matrix. We follow [57] and construct the triad (p, q, N) by a rotation 
around the vector e x with some constant inclination angle io, 

10 





cosi sini | | q | . (87) 
— sini cosi 

Figure [Tlshows shows a representation of what has been done. We clearly see that the vector p coincides 
with e x I3. Next, we express the radial separation r in the orbital plane (e x ,e y ) and perform the rotation 
Equation ([57]) to move to the observer's triad and calculate r and v, 



r = r (p cos 4> + q cos io sin 4> + N sin io sin <f) , 

v = P (j" cos — ^0 sin cj)j + q (r<j) cos io cos cf> + r cos io sin + N (r<j) sin io cos <p + r sin io sin . (89) 

This provides the orbital contributions to the field. To compute the spin contributions to the radiation field, 
we also expand the spins in the orbital triad, 

51 = Xi e z = Xi (Ncosio - q sini ) , (90) 

52 = X2e z = X2 (Ncosio - q sini ) . (91) 
We also need to know how h x and h + are extracted from the TT part. This is done via following projections: 

hx = ^(Piq j +P j q i )Vkl Tij h5 r i (92) 

h + = \ (PiPj - QjH) h l7 , (93) 



where V k , is the usual TT projector onto the line-of-sight vector N, 



1 



and we define 



= (SI - N l N k )(5f - WNi) - - N*W)(5 kl - N k N t ) , (94) 

V\f = ^{q iPj +p iqj ), (95) 

V^^liPiPj-HQj), (96) 
which are unaffected by the TT projection operator. 



§ In reference 1571 the caption for Figure 2 should be made precise. The plane of the sky meets the orbital plane at e x for 
T = only. Generally, at e x = p the plane of the sky meets the invariable plane. 
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Figure 1. The geometry of the binary. We have added the observer related frame (p, q, N) (in dashed and 
dotted lines) with N as the line-of-sight vector. N points from the origin of the invariable frame (e x , e y , e z ) 
to the observer. Note that the orbital angular momentum h lies on the e z axis and so do the spins. N is 
chosen to lie in the e y -e z -plane, and measures a constant angle io (associated with the rotation around e x ) 
from e z , such that p = e x , and this is the line where the orbital plane meets the plane of the sky. The angle 
io is also found between the vector q, itself positioned in (e y ,e z ), too, and e y . 



The above expressions, Equation (f86j) - (|96|) . enable us to compute all the considered contributions to the 
radiation field polarisations. Following [62] and [63], we list the lowest order contributions to the gravitational 
waveform in harmonic coordinates. These are the PP contributions to 2PN, including the NLO-SO and 
LO-S1S2 terms. We also add the terms emerging from the gauge transformation from ADM to harmonic 
coordinates. 



lTT 



2?7 



4T p + edT )pp 



2c t (l)SO 



2 f (l)PP 



4,(2) PP 



3 ^ S o^- 5)S ° + ^ 2 a SlS2 d- )SlS2 



e 2 6 a so ^ PP+S ° + ^ "so g d°- 5+1) PP+S ° + . 4 s d" +2) PP+PP 



(97) 



Those terms in the last line of the above equation, labeled "g", denote corrections coming from the gauge 
transformation from ADM to harmonic coordinates to the desired order |35t IM] . | Appendix B| gives deeper 
information about how velocities, distances and normal vectors change within this transformation. We find 
it convenient to give a hint to their origin by putting the GW multipole order and the order/type of the 
correction in the label, for example "(0 + 1)PP + SO" is the first Taylor correction of the "Newtonian" (PP) 
quadrupole moment where the coordinates are shifted by a 1PN SO term. 

According to Equation and Equation (|M|) one can define the projected components of the £ via 



(order) 



(order) 



_ •p( x )<r (order) 
ij Sij 

_ <p(+M° rdcr ) 
ij Sij 



where the "cross" and "plus" polarisations read 

>x,+ 



>(x,+) 



(98) 
(99) 

(100) 
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f (0.5) PP 
Sx.+ 



(1)PP 



x,+ 



Sm \ ,1 
— \ 3 N ■ n 12 - 

m r 



o-p( x >+) - f-p( x ^+) 



(N-v) 



:-p(x.+) _ o-p(x,+) 



1 



(1 - 3 V ) 



1 



(N-n^f-l ( 3v 2 -15r 
r 



,-.2 + 7 ]\ p(x.+) + 30rp(x,+) _ 14 p(x,+) 

f J 7171 nV VV 



(N-n 12 )(N-v) 



12r"P, m ' 327-' m , 



3(1 - 3?/)v 2 - 2(2 - 3?7) 



(x,+) 
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1 



(x,+) 



+ (N-v)' 



6"P, 



(x,+) 



-2-P, 



(x,+) 



p(x ) +) +4 I ) , (5 + 3 r?) p(x, + ) 



3(1 - 37?)f 2 - (10 + 3??)v 2 + 29- 



■p(> 



(i..5) pp _ Sm ( 1 



— — (1 - 27,) (N-n 12 r- 
m 12 r 



45v 2 - 105f 2 + 90- ] rV { ^ +) - 96r7>i*' +) 
r , 



42 v 2 - 21 Of 2 



- (N-n 12 ) 2 (N-v) - 

r 

-(N-n 12 )(N-v) 2 i 

r 



t ) -p(x,+) 



27v 2 - 135f 2 + 84M Vi*- +) + 336fV^< +) - 172Pi*'+) 



(N-v) 



1 



4_<p(x,+) _ 24T>( > 



^(N-n 12 ) - 
12 r 



(69 - 66??)v 2 - (15 - 90??)f 2 - (242 - 24??)- 



fp(x,+) 



(66 - 36??)v 2 + (138 + 84??)r 2 
1' 



(256 - 72??)- 



1 

12 



(N-v) 



p(x.+) + (i92 + 12r?)r^'+) 



(23 - 10??)v 2 - (9 - 18??)f 2 - (104 - 12t?)- 



p(x,+) 



+ 40t?) -fV^- 



(12-60t?)v 2 - (20 - 52??)- 



d 2) + P = ^(1 - 5t? + 57, 2 )| 2 40 (N • v) 4 7e + > - (N • n 12 ) 4 



■p(> 



90(v 2 ) 2 + (318^ - 1260f 2 )v 2 + 344^ + 1890f 4 



1 



-2310-f 2 m*' +) 



^1620v 2 + 3000^ - 3780f 2 Wi,*' +) - f 336v 2 - 1680f 2 + 688-^\V^' +) 



(N-n 12 ) 3 (N-v)i 

r 



^1440v 2 - 3360f 2 + 3600^ t^* ' +) 



1 



1608v 2 - 8040f 2 + 3 
120 (N-v) 3 (N-n 12 

^2 1 



>(x,+) 



- 3960r7> 



(x,+ 



(N-n 12 ) z (N-v) 



396v 2 - 1980f 2 + 1668- )V^- +) + 6480f"Pi*' +) 



36007> ( *< +) 



_ _L (N-v) 2 

30 v ' 



(87 - 315?? + 145t? 2 )v 2 - (135 - 465?; + 75t? 2 
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(101) 



(102) 



(103) 
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(289 - 9057/ + 1157/ 2 ) 



7?(x.+) 



240 - 660?/ - 240t/ 2 ] rVi y v < +) 



(30 - 270?/ + 630t/ 2 )v 2 - 60(1 - 67/ + 10t/ z 



>(x,+) 



-L(N-n 12 )(N-v)i 



(270- 11407/ + 1170t/ 2 )v 2 



(60 - 450?y + 900?? 2 )f 2 - (1270 - 3920t/ + 360t/ 2 )- 



;-p(x,+) 



(186 - 810t/ + 1450t? 2 )v 2 + (990 - 2910t/ - 930?7 2 )f 2 



(1242 - 4170t/ + 1930t/ 2 )- 

r 



+ 



1230 - 38107/ - 90t/ j 



-p(x,+) 

^■p(x,+) 



+ — (N-n 12 ) 2 i 
60 v ; r 



(117 - 4807/ + 540t/ 2 )(v 2 ) 2 - (630 - 2850?/ + 4050t? 2 )v 



(125- 740t/ + 900?/ 2 )iv 2 
r 



(105 - 1050?? + 31507/ 
(1048 - 31207/ + 240t/ 2 



2\ -4 



2-\ ^,2 



- (2715 - 8580t/ + 12607/ 2 )- r 

-p(x,+) 



(216 - 1380t/ + 4320t/ 2 )v 2 + (1260 - 3300?? - 36007/ 2 )f 2 



(3952 - 128607/ + 3660t/ 2 



(x,+) 



(12 - 1807/ + 1160r/ 2 )v 2 + (1260 - 3840?/ - 780r? 2 )f 2 



(664- 2360?y+ 1700t/ 2 ) 



' 11 



x,+ 



re 



1 

60 



r 

2w,,2\2 



(66- 15t/ - 1257/^)(v 2 ) 



+ (90 - 180r7 - 480?/ 2 )vV - (389 + 1030?/ - 110?/ 2 )- v 2 

r 

+ (45 - 225t/ + 225?/ 2 )t ; - 4 + (915 - 1440t/ + 720?/ 2 )- f 2 

r 



(1284 + 1090t/W 



_7?(x,+) 



(132 + 5407/ - 580?/ 2 )v 2 + (300 - 11407/ + 3007/ 2 )f 2 



(856 + 4007/+ 7007? 2 )- 



(45-315?/ + 585t7 2 )(v 2 



(270 - 307/ + 270?/ 2 )ir 2 
r 



i/.-p(x,+) 



(354 - 210t/ - 550t/ 2 ) - v 2 
r 
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,(l)SO 



(638 + 1400t? - 130?7 2 



•p(x>+) L 



{ [r^ x ' +) (A x N)* nij + y/T=4j[v£' +) (S x N)' n{ 2 ] } , 
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(104) 

(105) 



(1.5) SO _ 
x,+ — 



6P^' +) [v;A,nia]-6r 



^ x <+> (Axn 12 )«4 



+ 4 

+ 7? 

— 6r 



6p(^+)[v;S,n 12 ] 



(S x N) 1 n{ 2 (6r(N • n 12 ) - 4(N • v)) 



p| x ' +) (Sxn 12 )V L2 



+ 7/ U 



,(x,+) 



(S x n 12 ) 1 ^ - 4(N • n 12 ) 7>< x - +) (S x N) 1 v 



+ (277 + 4) pf x < +) (Sxv) 1 ^ 



(106) 



tf (2)SiS 2 3 ?7 , 

£x = TXiX2COs(7 )sm(2^) , 



t (2)SiS 2 



3?7 

4r 3 



X1X2 (cos (2io) cos(2</>) + 2 sin 2 (i ) + 3 cos(2^)) 



(107) 
(108) 



The remaining contributions are the gauge dependent terms. Explicitly, they read 



gf (0+l)PP+SO _ 
Sx,+ 



3V^ [v; S, n 12 ] + 2 V^ +) (S x v)* nj 2 + 2 ^ x ' +) (S x n 12 



(109) 



g ^(o.5+i)PP+so = Sm -q \ (x , +) 



2r 2 



-15r(N • n 12 )[v; S, n 12 ] - 3r[N; S, v] 



3(N-v)[v;S,n 12 ] - -[N;S,n 12 ] 



(N • n 12 ) 



— 6f 



P^ +) (Sxvyn{ 2 +187^+)[v;S,n 12 ] 



6 r 
r 



^• +) (S x n 12 ) J nj 2 + 6 P^' (S x v) 4 v 



,(x,+) 



+ (N • v) 



(S x v)* n{ 2 l + 4M x ' +) (S x m,) 1 



• 6P( X <+>[N; S, v] + 2^ X '+'[N; S, n 12 ] 



(110) 



g,(0+2)PP+PP _ 1 j-p(x,+) 
4. -p(x,+) 



1 /o-2 7 2N 2(577 - 1) 

-7? (3r - 7v ) 

77 (l7r 2 - 13v 2 ) + ™H±± 



4. -p(x,+) 



577 (v 2 - llf 2 ) I277+I 



4r 



2r 2 



(111) 



Motion and gravitational wave forms of eccentric spinning compact binaries 20 

Equation (jllip shows total agreement with the transformation term in Equation (A2) of [55]. The next block 
of equations evaluates the scalar products of vectors and projectors containing the spins. First, we list those 
with the total spin S = Si + S2. For those terms with A = Si — S2 instead of S, simply replace (S — > A) 
on the left hand side and (xi + X2) (xi — X2) on the right. The used abbreviations are given by 



[S;ni 2 ,v] =Mxi+X2), (112) 

[N;S,m 2 ] =(xi+xa)sin(io)cos(0), (113) 

[N;S,v] = (xi+X2)sin( io )(fcos(0)-0rsin(0)), (114) 

7>4 X V(S X n 12 )' = ^(xi+Xa) cos (i ) {rcos(20) - <^ sin(20)} , (115) 
^'(S x ma) 1 = 1 (xi + Xa) |-^r (cos (2i ) + 3) cos(20) + 2<£r sin 2 (*„) 

-r(cos(2t ) + 3)sin(20)| , (116) 

V\f n\ 2 (S x n 12 ) j = 1 (xi + X2) cos (i ) cos(20) , (117) 

^ < 2 (S x ni 2 ) J - - ^ (Xi + Xa) {cos (2i ) + 3} sin(20) , (118) 

7^ x V(Sxv) J = - -(xi + xa) cos (to) {2^rrsin(20) + oos(20)(0r- r)(</>r + r)} , (119) 

pf'n'iSxv) 1 =-(xi+X2)(cos(2i o ) + 3){sin(20)(^r-r)(^r + f)-2^rrcos(20)} , (120) 

pJ'nUSxvy = i(xi+X 2 )cos(i o ){fcos(20)- 0rsin(20)} , (121) 

nl 2 (S x v) J = g (xi + X2) {-0r cos(20) (3 + cos(2i )) - 20r sin 2 (i ) 

-rsin(2^)(3 + cos(2i ))}, (122) 

(S x N) V 12 = - l(xi + xa) sm(t ) cos(0) , (123) 

7^ x) (S x N)V 12 = - J(xi + xa) sin(2* ) sin(^) , (124) 

V\p (S x N) V = - I ( Xl + xa) sin(z ) (f cos(^) - & sin(0)) , (125) 

P§°(SxN)V = -i(xi+X2)sin(2i o )(0rcos((/))+rsin(0)). (126) 
The spin-independent projections and the ratio of the difference to the sum of the masses read 

(N • n 12 ) = N z n\ 2 = sin(io) sin(<£) , (127) 

(N-v) = N.y =sin(t o )|r0cos(^) +fsin(0)| , (128) 

v 2 = v*Vi = r 2 <j) 2 + r 2 , (129) 

^ = "P^nUnit = cos (to) sin(0) cos(0) , (130) 

V™ =r i ^ ) v i v j =icos(i ){sin(2(/))(f 2 -(/) 2 r 2 ) +20frcos(2(/))} , (131) 

^ =^ X) <2^' =^cos(io){^rcos(2^+rsin(2^)} , (132) 

= ^^{2 = \ {cos 2 (0) - cos 2 (to) sin 2 (0)} , (133) 

Vit ] = V^ ) v i v j = - {(rcos(0) - <^rsin(0)) 2 - cos 2 (t ) (</>rcos(0) + r sin(c/>)) 2 } , (134) 

= P^nja^ = g (cos (2i ) + 3) sin(20) - 4f cos 2 (i ) sin 2 (</>) + 4f cos 2 (c/>)} , (135) 
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mi — 777,2 



(77 



(136) 



5m 

777 

In the expression for the emitted gravitational wave amplitudes, Equation ([57)) , R' is the rescaled distance 
from the observer to the binary system, 

G?77 



R 1 = R 



(137) 



We note that it is very important that R' has got the same scaling as r in order to remove the physical 
dimensions. The common factor c~ 4 of hj^ will be split in c~ 2 for the distance R 1 and c~ 2 for the , in 
order to make all terms dimensionless. Also note that the Equations (|100[) - (|135p in our special coordinates 
are valid only when h is constant in time. In the non-aligned case, additional angular velocity contributions 
kick in and the expressions become rather impractical. From reference [B], the reader can extract explicit 
higher-order spin corrections to the Newtonian quadrupolar field for the case of quasi-circular orbits. 



7.2. Dynamical orbital variables as implicit functions of time 

We are now in the position to compute the time domain gravitational wave polarisations with the help of our 
orbital elements, to be expressed in terms of conserved quantities and the mean anomaly, which is an implicit 
function of time. Using Equations (|70|) . (|71|) . (fT2|) and (fT3"f one can express the quantities r, f, (f>, <f> (which 
are used in the radiation formulas) in terms of the eccentric anomaly 77, other orbital elements and several 
formal 2PN accurate functions. The most compact quantity is r which is given by Equation (|70[) . namely 



r(u) = a r (l — e r cosu) 
(1 — et cos u) 



2\E\ 



\E\ 



\E\ 



2(1 — et cosu) 
1 



{(1 — e t cos7i)(9 — 577) + 677 — 16} 



(1 - e t cos u) 1 4 (ef — 1) 



+ 277 (-377 (e 2 t - 1) + Me 2 - 56) + 68 

|£|3/2 



(1 - e t cos7t)(?7(7?7 - 58) + 1) (e 2 - l) 
6((1 - e t cos u) - l)(2?7-5) 



0" 



e Sa s 



(1 — e t cosu)-\/l — e 



{2WI-477 (xi - xa) - V5(V - 2) (Xi + X2)} 



(138) 



+ 



(1 — et cos 77) 

(Xi -X2)' 



+ (Xi + X2) ! 



Oh' ((Ai + A 2 )(2t7 - 1) - (Ai - \ 2 )VT^) na SlS2 
2(ef-l) e?-l 

a s 2 ((Ai +A 2 )(2?7- 1) - (Ai - Az)^! - 4^) 77a SlS; 



2(e?-l) 



e?-l 



a s 2 (xi ~ X2) (xi + X2) 



^|£;| 5 / 2 



(Ai - A 2 )(27 ? - 1) - (Ai + \ 2 )^/l^ 



2\/2(l — e f cos 77) 

J 16((1 — e t cos 77) — 1) 



(e? ~ 1) 



((77 - 8)77 + 6) (xi + xa) - 2(77 - 3)v/1-4t7( Xi - Xa ) 



+ 



3(77(277 - 15) + 12) (xi + X2) - (1377 - 36)v/1-4t7 ( X i - X2) 



(1977 " 42) (Xi - X2) - (77(13/7 - 50) + 42) ( X i + X2) 



(139) 
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Using expression (|138|) . we calculate the derivative via chain rule, given by 
. . . dr du 

r(u) = — = na r e r sin it x 

y ' dudt 
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1 _ e %{ e <t> ~ cosw) 

1 - et cos u + T v J —r. ^3 V Fv-u 

(1 — e^cosu)^ 



1-ei 



1 — cos u 



- 1 



(140) 



e t y / 2|£|sini 
1 — et cos u 



1 + e 2 \E\ -(1 - 3r?) + e 4 \E\ 2 —(23 + 77(4777 - 98)) 



32 



6(2r?-5)Vl 



1(77 - 20)77 - 30 77(77 + 4) (e 2 - l) 



(e 2 — 1) (1 — e t cosm) (1 — e t cosw) 2 2 (1 — e t cos u) 3 
+ a so <5 x 

• 4^2 (-2^1— 4^(77 - 3) (xi - X2) + (fa - 8)77 + 6) (xi + X2)) 



e 4 |£ 



5/2 



(e 2 — 1) (1 — e t cosu) 
(yT^(17T7 - 40) (xi - X2) + ((51 - 877)77 - 40) (xi + X2)) 

\/2 - 2e t 2 (1 - e t cost/,) 2 
y/1 - e 2 (VT^fa + 8) (xi - X2) + (8 - 1377) (xi + X2)) 
-\/2(l — e t cosu) 3 



(141) 



The final expression for <f> in terms of u is rather complicated. It is convenient to give a short expression and 
a description how to obtain it. From Equations (|72p and (|73p one can eliminate v to obtain 



(7JI = (p 





ll + e 4 


4. 


— < 2 arctan 




■ tan — 


2tt\ 


.V 1 - e ^ 


2 



2^1 - e| sin(u)(e0 - cos(tt)) 



(e^ cos(u) - l) 2 



.^1 - e 2 sin(u) ((e 2 - 4) cos(2?i) - 7e 2 + 12e cos(u) 
2(e0 cos(m) — l) 3 



(142) 



Using the chain rule once more one gets an expression for the angular velocity via 4>{u) 
(d(j)/dv) (dv/du) (du/dt), symbolically 



4>{u) 



1-e 



P (1 — cosu) 

^ - 4e cos(u) - (e^ - 2) cos(2tj)) 



1 + Q2v 



G3v 



(e^ cos(Tt) - l) 2 

(30e 3 - 45e 2 cos(ii) - 18 (e 2 - 2) e cos(27t) + 3 ^3e 2 - 4) cos(3u)) 



4(e0 cos(u) — l) 1 



1 — et cos it + J- v 



1 - e^(e0 - cosw) 
(1 — cos u) 2 



J v—u 



I -el 



1 — CA, COS 7i 



(143) 



Again, P can easily be computed with the help of the already known definition n = 2ir/P and Equation 



8. Conclusions 



In this paper we presented a quasi-Keplerian parameterisation for compact binaries with spin and arbitrary 
mass ratio. We assumed that the spins are aligned or anti-aligned with the orbital angular momentum and 
restricted ourselves to the leading-order spin-spin and next-to- leading order spin-orbit, as well as to 2PN 
point particle contributions. 
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The conservation of alignment for all times holds if the alignment is assumed at the initial instant of 
time. It turned out that the effects of the spins do not destroy the polynomial structure of the integrals for 
both the angular and the radial variables, for which the standard routine is valid [23l [59l [65] and enabled us 
to give a fully analytic prescription for the orbital elements in terms of the binding energy, the mass ratio 
and the magnitudes of the angular momenta. 

Furthermore we provided PN extended formulae for the radiation polarisations in analytic form derived 
from the results of [SH HH] due to the currently highest available order in spin. The missing SO part at 2PN 
will be given in a forthcoming publication. 

An outstanding question is the stability of the spin configurations. If we assume that the spins have 
tiny differences in their directions, it is interesting to know if the enclosed angles will grow secularly or will 
oscillate in an unknown manner. This will be task of a further investigation, as well as the inclusion of 
additional higher order spin Hamiltonians. 

Another task to be tackled is the radiation reaction. It is possible to include the conservative 
contributions of the spin to the orbital motion into the equations of the far-zone energy and angular 
momentum flux expressions. The goal is an equation of motion for the orbital elements to be obtained 
in an adiabatic approach. 

Acknowledgment s 

We thank Jan Steinhoff and Steven Hergt for many useful discussions. This work is partly funded by the 
Deutsche Forschungsgemeinschaft (DFG) through SFB /TR7 "Gravitationswellenastronomie" and the 
Research Training School GRK 1523 "Quanten- und Gravitationsfelder" and by the Deutsches Zentrum fur 
Luft- und Raumfahrt (DLR) through "LISA Germany". 



[1] S. Rowan and J. Hough, "Gravitational wave detection by interferometry (ground and space)," Living Rev. Relativity 3 

(2000) 3. |http : //www . liv ingreviews . org/lrr-20 00-3| 
[2] G. Faye, L. Blanchet, and A. Buonanno, "Higher-order spin effects in the dynamics of compact binaries. I. Equations of 

motion," \Phys. Rev. D 74 (2006) 104033[ |arXiv:gr-qc/0605l39| 
[3] L. Blanchet, A. Buonanno, and G. Faye, "Higher-order spin effects in the dynamics of compact binaries. II. Radiation 

field," \Phys. Rev. D 74 (2006) 104034[ |arXiv:gr-qc/0605i40) 
[4] L. Blanchet, A. Buonanno, and G. Faye, "Erratum: Higher-order spin effects in the dynamics of compact binaries. II. 

Radiation field," \Phys. Rev. D 75 (2007) 049903(E)) 
[5] L. Blanchet, A. Buonanno, and G. Faye, "Erratum: Higher-order spin effects in the dynamics of compact binaries. II. 

Radiation field," \Phys. Rev. D 81 (2010) 089901(E)) 
[6] K. G. Arun, A. Buonanno, G. Faye, and E. Ochsner, "Higher-order spin effects in the amplitude and phase of 

gravitational waveforms emitted by inspiraling compact binaries: Ready-to-use gravitational waveforms," 

\PTys. Rev. D 79 (2009) 104023||arXiv: 0810 .5336 [gr-qc]") 
[7] M. Hannam, S. Husa, B. Briigmann, and A. Gopakumar, "Comparison between numerical-relativity and post-Newtonian 

waveforms from spinning binaries: The orbital hang-up case," \Phys. Rev. D 78 (2008) 104 007[ 

[arXiv: 0712. 3787 [gr-q"cT| 

[8] P. Jaranowski and G. Schafer, "Radiative 3.5 post-Newtonian ADM Hamiltonian for many-body point-mass systems," 

\Phys. Rev. D 55 (1997) 4712-4722| 
[9] L. Blanchet, G. Faye, B. R. Iyer, and B. Joguet, " Gravitational- wave inspiral of compact binary systems to 7/2 
post-Newtonian order," \Phys. Rev. D 65 (2002) 061501 (R)) |arXiv:gr-qc/0105099) 
[10] L. Blanchet, G. Faye, B. R. Iyer, and B. Joguet, "Erratum: Gravitational-wave inspiral of compact binary systems to 7/2 
post-Newtonian order," \Phys. Rev. D 71 (2005) 129902(E7 



[11] P. Ajith, S. Babak, Y. Chen, M. Hewitson, B. Krishnan, A. M. Sintes, J. T. Whelan, B. Briigmann, P. Diener, 
N. Dorband, J. Gonzalez, M. Hannam, S. Husa, D. Pollney, L. Rezzolla, L. Santamaria, U. Sperhake, and 
J. Thornburg, "Template bank for gravitational waveforms from coalescing binary black holes: Nonspinning binaries 
\Phys. Rev. D 77 (2008) 104017||arXiv: 0710 .2335 [gr-qc] 



[12] P. Ajith, S. Babak, Y. Chen, M. Hewitson, B. Krishnan, A. M. Sintes, J. T. Whelan, B. Briigmann, P. Diener, 
N. Dorband, J. Gonzalez, M. Hannam, S. Husa, D. Pollney, L. Rezzolla, L. Santamaria, U. Sperhake, and 
J. Thornburg, "Erratum: Template bank for gravitational waveforms from coalescing binary black holes: Nonspinning 
binaries," \Phys. Rev. D 79 (2009) 129901(E)) 

[13] M. Mathisson, "Neue Mechanik materieller Systeme," Acta Phy s. Pol. 6 (1937) 163-200. 

[14] A. Papapetrou, "Spinning test-particles in general relativity. I," |Proc. R. Soc. A 209 (1951) 248-258] 

[15] E. Corinaldesi and A. Papapetrou, "Spinning test-particles in general relativity. II," |Proc. R.lS oc. A 209 (1951) 259-268 



Motion and gravitational wave forms of eccentric spinning compact binaries 



24 



B. M. Barker and R. F. O'Connell, "Derivation of the equations of motion of a gyroscope from the quantum theory of 

gravitation," \Phys. Rev. D 2 (1970) 1428-1435| 
B. M. Barker and R. F. O'Connell, "Gravitational two-body problem with arbitrary masses, spins, and quadrupole 

moments," \Phys. Rev. D 12 (1975) 329-335| 
P. D. D'Eath, "Interaction of two black holes in the slow-motion limit," \Phys. Rev. D 12 (1975) 2183-2199| 

B. M. Barker and R. F. O'Connell, "The gravitational interaction: Spin, rotation, and quantum effects — a review," 
| Gen. Relativ. Gravit. 11 (1979) 149-175| 

T. A. Apostolatos, "Search templates for gravitational waves from precessing, inspiraling binaries," 

\Phys. Rev. D 52 (1995) 605-620] 
L. E. Kidder, C. M. Will, and A. G. Wiseman, "Spin effects in the inspiral of coalescing compact binaries," 

\Phys. Rev. D 47 (1993) R4183-R4187||arXiv:gr-qc/9211025| 
T. Damour and N. Deruelle, "General relativistic celestial mechanics of binary systems. I. the post-Newtonian motion.," 
Ann. Inst. H. Poincare A 43 (1985) 107-132. 

G. Schafer and N. Wex, "Second post-Newtonian motion of compact binaries," \Phys. Lett. A 174 (1993) 196— 205| 

G. Schafer and N. Wex, "Erratum: Second post-Newtonian motion of compact binaries," \Phys. Lett. A 177 (19 93) 461(E) 
R.-M. Memmesheimer, A. Gopakumar, and G. Schafer, "Third post-Newtonian accurate generalized quasi-Keplerian 
parametrization for compact binaries in eccentric orbits," \Phys. Rev. D 70 (2004) 104011 1 [arXiv : gr-qc/0407049 

C. Reisswig, S. Husa, L. Rezzolla, E. N. Dorband, D. Pollney, and J. Seiler, "Gravitational- wave detectability of 
equal-mass black-hole binaries with aligned spins," Phys. Rev. D 80 (2009) 124026 ] [arXiv : 0907 ■ 04 62 [gr-qc] 

T. Bogdanovic, C. S. Reynolds, and M. C. Miller, "Alignment of the spins of supermassive black holes prior to 

coalescence," \ApJ 661 (2007) L147-L150| |arXiv : gr-qc/0703054l 
J. Steinhoff, S. Hergt, and G. Schafer, "Next-to-leading order gravitational spin(l)-spin(2) dynamics in Hamiltonian 

form," \Phys. Rev. D 77 (2008) 081501(R)] |arXiv: 0712. 1716 [gr-qc]"| 
S. Hergt, J. Steinhoff, and G. Schafer, "The reduced Hamiltonian for next-to-leading-order spin-squared dynamics of 

general compact binaries," \Class. Quant. Gray. 27 (2010) 1 35007 , arX iv: 1002 .2093 [gr-qc] | 
J. M. Martin-Garci'a, "xPerm: fast index canonicalization for tensor computer algebra," 

\Comp. Phys. Commun. 179 (2008) 597-603] |arXiv: 0803. 0862 [cs.SCJI 

J. M. Martfn-Garcia, xAct: Efficient Tensor Computer Algebra. |http : //www . xact ■ es/ 1 

P. Jaranowski and G. Schafer, "Third post-Newtonian higher order ADM Hamilton dynamics for two-body point-mass 
systems," \Phys. Rev. D 57 (1998) 7274-7291] |arXiv : gr-qc/9712075] 



Damour, P. Jaranowski, and G. Schafer, "Dimensional regularization of the gravitational interaction of point masses," 
\Phys. Lett. B 513 (2001) 147-15*5] |arXlv:gr-qc/010 5038 
E. Poisson, "Gravitational waves from inspiraling compact binaries: The quadrupole-moment term," 

\Phys. Rev. D 57 (1998) 5287-5290||arXlv:gr-qc/9709032| 
T. Damour, P. Jaranowski, and G. Schafer, "Hamiltonian of two spinning compact bodies with next-to-leading order 

gravitational spin-orbit coupling," \Phys. Rev. D 77 (2008) 64032 , ar Xiv: 0711 . 1048 [gr^qcT] 
J. Steinhoff, G. Schafer, and S. Hergt, "ADM canonical formalism for gravitating spinning objects," 

\Phys. Rev. D 77 (2008) 104018||arXlv: 0805 .3136 [gr-qc] | 
M. Levi, "Next to leading order gravitational spinl-spin2 coupling with Kaluza-Klein reduction," 

[arXiv: 0802. 1508 [gr-qc]"| 
K. S. Thorne and J. B. Hartle, "Laws of motion and precession for black holes and other bodies," 

\Phys. Rev. D 31 (1985) 1815-1837| 
S. Hergt and G. Schafer, "Higher-order-in-spin interaction Hamiltonians for binary black holes from source terms of Kerr 

geometry in approximate ADM coordinates," \Phys. Rev. D 77 (20 08) 104001] |arXiv: 0712 . 1515 [gr-qc] | 
S. Hergt and G. Schafer, "Higher-order-in-spin interaction Hamiltonians for binary black holes from Poincare invariance," 

\Phys. Rev. D 78 (2008) 124004||arXiv: 0809 .2208 [gr-qc]~| 
J. Steinhoff and G. Schafer , "Canon ical formulation of self-gravitating spinning-object systems," 

\Europhys. Lett. 87 (2009) 50004||arXiv: 0907. 1967 [gr-qc]"| 
J. Steinhoff and H. Wang, "Canonical formulation of gravitating spinning objects at 3.5 post-Newtonian order," 

\Phys. Rev. D 81 (2010) 024022| [arXiv: 0910 ■ 1008 [gr-qc]~| 
W. M. Tulczyjew, "Motion of multipole particles in general relativity theory," Acta Phys. Pol. 18 (1959) 393-409. 
W. G. Dixon, "Extended bodies in general relativity: Their description and motion," in Proceedings of the International 

School of Physics Enrico Fermi LXVII, J. Ehlers, ed., pp. 156-219. North Holland, Amsterdam, 1979. 
J. Steinhoff and D. Puetzfeld, "Multipolar equations of motion for extended test bodies in general relativity," 

\Phys. Rev. D 81 (2010) 044019||arXlv: 0909 .3756 [gr-qc] | 
W. G. Laarakkers and E. Poisson, "Quadrupole moments of rotating neutron stars," \Astrophys. J. 512~ (1999) 282-287 
|arXlv : gr-qc/9709033l 

J. Levin, "Gravity waves, chaos, and spinning compact binaries," \Phys. Rev. Lett. 8 4 (2000) 3515-3518 
|arXiv:gr-qc/9910040] 

C. Sohr, "Chaos in gravitativ-gebundenen Binarsystemen mit Spin.," diploma thesis, Friedrich-Schiller-Universitat Jena, 
2009. unpublished. 

P. A. M. Dirac, Lectures on Quantum Mechanics. Yeshiva University Press, New York, 1964. 
H. Goldstein, Classical Mechanics. Addison- Wiley Publishing Company, 2nd ed., 1981. 
P. Colwell, Solving Kepler's equation over three centuries. Willman-Bell, Inc., Richmond, VA 23235, 1993. 
B. Mikoczi, M. Vasuth, and L. A. Gergely, "Self-interaction spin effects in inspiralling compact binaries," 
\Phys. Rev. D 71 (2005) 124043||arXlv:astro-ph/0504'538] 
[53] J. Majar and M. Vasuth, "Gravitational waveforms from a Lense-Thirring system," \Phys. Rev. D 74 (2006) 124007 
|arXiv:gr-qc/0611105| 



Motion and gravitational wave forms of eccentric spinning compact binaries 



25 



[54 
[55 
[56 
[57 
[58 
[59 
[60 
[61 
[62 
[63 

[64 

[65 

[66 



M. Vasuth and J. Majar, "Gravitational waveforms for finite mass binaries," \Int. J. Mod. Phys. A 2 2 (2007) 2405 -2414] 
|arXiv: 0705. 3481 [gr-qc] | ^ 

J. Majar and M. Vasuth, "Gravitational waveforms for spinning compact binaries," \Phys. Rev. D 77 (2 008) 104005 
|arXiv = 0806.2273 [gr-qc] | 

L. A. Gergely, "Spinning compact binary inspiral: Independent variables and dynamically preserved spin configurations," 
larXiv:0912.0459l 

M. Tessmer, "Gravitational waveforms from unequal-mass binaries with arbitrary spins under leading order spin-orbit 

coupling," \Phys. Rev. D 80 (2009) 124034| |arXiv: 0910. 5931 [gr-qc] | 
Z. Keresztes, B. Mikoczi, and L. A. Gergely, "The Kepler equation for inspiralling compact binaries," 

\Phys. Rev. D 72 (2005) 104022] |arXiv:astro-ph/0510602| 
T. Damour, A. Gopakumar, and B. R. Iyer, "Phasing of gravitational waves from inspiralling eccentric binaries," 

\Phys. Rev. D 70 (2004)1) 64028 ar Xlv:gr-qc/0404128| 
C. Konigsdorffer and A. Gopakumar, "Phasing of gravitational waves from inspiralling eccentric binaries at the 

third-and-a-half post-Newtonian order," \Phys. Rev. D 73 (2006) 124012| |arXiv : gr-qc/0603056] 
L. A. Gergely, Z. I. Perjes, and M. Vasuth, "Spin effects in gravitational radiation backreaction. Ill: Compact binaries 

with two spinning components," \Phys. Rev. D 58 (1998) 12400T] |arXlv: gr-qc/9808063 
L. E. Kidder, "Coalescing binary systems of compact objects to (post)'" 3 ^-Newtonian order. V. Spin effects," 

\Phys. Rev. D 52 (1995) 821-847] |arXiv:gr-qc/9506022] 
A. Gopakumar and B. R. Iyer, "Gravitational waves from inspiraling compact binaries: Angular momentum flux, 

evolution of the orbital elements, and the waveform to the second post-newtonian order," 

\Phys. Rev. D 56 (1997) 7708-7731||arXiv:gr-qc/01l6Too] 
T. Damour, P. Jaranowski, and G. Schafer, "Equivalence between the ADM-Hamiltonian and the harmonic-coordinates 

approaches to the third post-Newtonian dynamics of compact binaries," Phys. Rev. D 63 (2001) 044021 

arXiv : gr-qc/0010040 

C. Konigsdorffer and A. Gopakumar, "Post-Newtonian accurate parametric solution to the dynamics of spinning compact 
binaries in eccentric orbits: The leading order spin- orbit interaction," \Phys. Rev. D 71 (2005) 024039| 
arXiv : gr-qc/0501011 

T. Damour and G. Schafer, "Lagrangians for n point masses at the second post-Newtonian approximation of general 
relativity," \Gen. Relativ. Gravit. 17 (1985) 879-905| 



Appendix A. Integrals 

For the sake of completeness we give the results of the definite integrals /„ and I' n for different n: 

A = -^=> (A.2) 

I' 2 = 2tt, (A.3) 
=*(*_+*+), (A.4) 



V A = -7T (3*L + 2s_,s+ + 34) , (A.5) 

I' 5 = £7r(s_ + s+) (5s 2 _ - 2s-s+ + 54) . (A.6) 
The more complicated integrals with boundary s in terms of u, v, e r and a r are given by 



lo = a 2 r \J\ — e 2 ;. (u — sinu), (A. 7) 

h = a rV /l - eju, (A.8) 

h = v, (A.9) 

/. - ' + f: Sin Jf ; (A.10) 
a r (1 - e$) 

2(2 + e^)w + 8e r sinu + e 2 S in(2w) 

h AaUl-eD* ' (A - 11} 

6 (2 + iel) v + 9e r (4 + el) sin v + %el sin <2v) + el sin (3v) 

h = " 12a3(l-e?)3 • (A ' 12) 



Motion and gravitational wave forms of eccentric spinning compact binaries 



2G 



Appendix B. Coordinate transformation from ADM to harmonic 

From section IV of |35j and from |64| . we collect the contributions for the coordinate transformation from 
ADM to harmonic coordinates for spinning compact binaries, including LO effects of spin-orbit interaction 
and 2PN PP contributions. Let Y a label the harmonic position of the a-th particle as a function of the ADM 
positions x b , momenta p b and spins S b . Then, to 2PN order, the transformation reads in their notation 



Y a (x b , p 6 ) = x a + e 2 YS°( Xb , p b , S b ) + e 4 Y 2PN (xfc, p 6 ) 
with 



Yf°(x h , Ph ,S b ) 



S a x p a 



2ml 



Y 2PN (x a ,p a ) =Gm 2 



8 m\ 8 m\ 



(B.l) 
(B.2) 



5 p 2 1 (ni2 • P2) Gmi ( 7 1 m 2 



7*12 V4 4 mi 



ni2 



1 (n 12 ■ p 2 ) Pi 
2 



m 2 



in 1 



7 (n 12 ■ p 2 ) p 2 

4 



m 2 



m 2 



(B.3) 



where Y 2PN (x a , p a ) is simply obtained by exchanging the particle indices (1 -f-> 2). We find it very important 
to mention some of the rules to obtain the relative separation vector with the scaling introduced in this 
paper. The above equations are not given in relative coordinates. Thus, we scale every S a with m 2 . Next, 
we subtract Yi from Y 2 , setting p 2 = pi = — p for the centre-of-mass frame and scale p with /1 as in 
Equation (J8J) to get a dimensionless momentum. Finally, we divide the obtained separation vector with G m 
and obtain the separation in terms of the linear momentum and the ADM spin momenta. 

The harmonic velocity is obtained just by plugging the harmonic positions in the Poisson brackets with 
the total Hamiltonian, 

v harm = ^harm^ADMj ^ ^ ^ 

see [SS] for the effects of the 1PN time shift part of the coordinate gauge transformation to the velocities 
and positions. The linear momentum p can then be expressed in terms of the velocity perturbatively. It is 
important to express p in terms of the ADM velocity first and then to plug it into the expression for v harm 



afterwards. To 2PN order, the radial separation, the velocity and the unit normal vector, r 
harm transform due to 



v 



ni2 

„harm 



x + ^e 2 5i] (S x v) + e 
v-e 2 ^(Sx n 12 ) 



127?+ 1. 

4r 



1 ( 
-n i2 - -77 (ni2 



5v^ 



18rv 



e 4 j|- (fn 12 (3f 2 - 7v 2 ) + v (l7r 2 - 13v 2 )) + ((21r? + l)v - (19r? + 2)rn 12 )| 



„harm 



-e 2 5r] [S;n 12 ,v] + e 4 



-77 (5v 2 - 19f 2 ) 



3?7 



ni2 



ha nil 



= n 12 



2r 



{n i2 [S;ni2,v] + (S x v)} 



6 -^r{ rn ^-v} , 



where every quantity on the right hand side is written in ADM coordinates. Note that 

ni2 [S; ni 2 , v] + (S x v) = (1 - n i2 (8 m 2 ) (S x v) 
is the part of (S x v) which is orthogonal to ni 2 . 
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